ABSTRACT 

'X  transfer  function  approach  is  applied  to  the  problem  of 
lightning  coupling  into  internal  aircraft  avionic  circuits.  This 
approach  allows  a systematic  procedure  for  determining  transient 
levels  Induced  on  interior  aircraft  circuits  due  to  exterior  inter- 
ference fields  such  as  those  caused  by  nearby  or  direct-strike 
lightning.  Unique  in  the  approach  used  here  is  the  use  of  a tri  * 
angular  patch  model  for  the  exterior  surface  of  the  aircraft  with 
all  points  of  entry  (PO£s)  closed,  ,^A  moment  method  procedure  is 
then  used  to  find  the '■"short  circuit**  skin  current  in  the  frequency 
domain.  This  current  interacts  with  all  the  POEs,  described  by 
either  measured  or  theoretically  calculated  transfer  functions, 
in  a secondary  interaction  problem  to  create  interior  penetrating 
fields.  These  interior  fields  are  used  to  derive  distributed 
voltage  and  current  sources  which  excite  a transmission  line 
model  of  a wire  harness  or  cable  bundle.  Further  levels  of 
coupling  are  indicated  along  with  their  concomitant  transfer 

functions. 
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INTRODUCTION 

THE  GENERAL  PROBLEM  of  scslytically  deter- 
mining the  susceptibility  of  interior  air- 
craft avionics  circuitry  to  upset  or  damage 
caused  by  an  external  electromagnetic  threat 
is  Inherently  complicated.  To  facilitate 
the  analysis,  two  principle  assumptions 
concerning  the  electromagnetic  coupling 
phenomenon  are  made  at  the  outset: 

* Linearity  of  material  media 

* Interaction  level  decoupling 

The  meaning  of  the  first  assumption  is 
obvious  and  test  results  show  that  it  is 
generally  valid  for  bulk  cable  current 
responses.  The  second  assumption  means 
that  the  aircraft  has  no  effect  on  the 
source,  the  interior  compartments  have  no 
effect  on  the  external  skin  currents,  and 
the  interior  circuits  have  no  effect  on  the 
interior  fields.  It  is  easy  to  construct 
pathological  cases  for  which  the  second 
assumption  is  theoretically  invalid. 
Measured  results,  however,  indicate  that 
for  a great  many  cases,  it  is  a reasonable 
assumption. 

Figure  1 illustrates  a matrix  of 
possible  interaction  levels  where,  at  each 
level,  a source  Interacts  with  a 
characteristic  to  produce  a result.  Further 
simplifying  assumptions  concerning  the 
threat,  the  external  physical  detail  of  the 
aircraft,  and  the  interior  detail  of  the 
circuit  in  queatlon  must  still  be  made  in 
order  to  make  the  problem  tractable.  Errors 
introduced  at  each  characteristic  model 
level  combine  to  csst  accumulative 
uncertainty  at  each  succeaaive  result 
level.  Hopefully,  careful  comparison 
between  computations  and  measured  results 
will  aid  in  refining  characteristic  models 
and  reducing  result  error. 

Using  the  first  two  assumptions,  *>nd 
working  in  the  frequency  domain,  one  my 
express  the  results  of  level  D in  Fig.l  by 
the  following  equations: 
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where : 

Vg(f),  Ig(f)  " Voltage  and  current  trans- 

fer function  at  input  of 
susceptible  device  (avionics 

box)  f 

D(f)  ■ Threat  driving  function 

spectrum 

TA(,(f)  » Band-limited  approximation 

to  Impulse  response  of 
aircraft 

TpOE^)  " Polnt-°f-entry  (POE)  trans- 

transfer function 
Z - Material  shielding 
transfer  impedance 
Yj  * Joint  admittance 

T^(f)  ■ Transfer  functions  account- 
ing for  compartment  detail  I 


_ V I 

Tc  ’ Tc 

l 

T_(f) 

- Transmission  line 

coupling 

1 

transfer  function 

including 

cable  shielding 

1 

Vf> 

” Avionics  box 

protection 

function  including  penetra- 
tion and  circuit  protection 
functions 

This  juxtaposition  of  transfer 
functions  is  illustrated  plctorlally  in 
Fig.  2.  Further  detail  concerning  the 
above  transfer  function  approach  to 
electromagnetic  coupling  applied  to 
aircraft  is  given  in  [1]*  and  (2]. 
Examples  of  some  of  the  above 
characteristic  transfer  functions  are 
illustrated  plctorlally  in  Fig.  .. 

DESCRIPTION  OF  THE  METHOD 

The  double  exponential  time  function 

D(t)  - A(e'°t  - e'Bt)  (2) 

is  often  used  to  model  lightning  waveforms 
[3]  vhere  the  psrameters  a snd  0 may  be 
adjusted  to  change  the  rise  time  1/0  and 
fall  time  1/a  of  the  lightning  return 
stroke.  This  simple  double  exponential 
model  is  assumed  to  sufficiently 
characterise  the  behavior  of  the  current  in 

•Numbers  in  brackets  Indicate  references  at 
end  of  paper. 
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the  lightning  channel.  In  this  model,  the 
lightning  attachment  to  the  aircraft  is  well 
underway  and  the  source  of  the  channel  is 
treated  as  a very  high  impedance  current 
generator  so  that  the  aircraft  has  no 
loading  effect  on  the  channel. 

When  a conducting  aircraft  becomes 
part  of  the  path  of  a lightning  channel 
current,  a large  charge  displacement  takes 
place  at  certain  "junction"  points  on  the 
aircraft  where  the  channel  enters  and 
exits.  These  points  are  deemed  "entry"  and 
"exit"  points,  where  the  conventional 
reference  direction  of  current  is  implied 
to  denote  charge  transfer.  Here  we  assume 
that  the  entry  and  exit  points  are  given 
and  that  they  coincide  with  nodes  in  the 
triangular  surface  patch  model  {4,  S,  and  6] 
of  the  aircraft  skin.  This  will,  in  fact, 
be  the  case  when  modelling  an  experimental 
setup  where  current  injection  points  are 
well  defined  [7].  Further  details 
concerning  the  computation  of  direct  strike 
skin  currents  may  be  found  in  [6]  and 
details  concerning  direct  strike  test 
methods  may  be  found  in  [8,  9,  and  10]  for 
operational  aircraft. 

A widely  used  and  readily  available 
data  base  which  models  aircraft  skin  by 
surface  patches  is  the  so-called  NASTRAN 
data  base  [11] . An  example  of  this  is 
shown  in  Fig.  3 for  the  Black  Hawk 
helicopter.  In  order  to  solve  for  the  skin 
currents  >/  the  method  of  moments,  a reduced 
data  bare  consisting  of  triangular  surface 
patches  is  created.  Examples  of  these  are 
shown  i.n  Figs.  4 through  7 for  the  Black 
Hawk  helicopter  and  F-14  fighter  aircraft. 
Here  we  have  picked  various  direct-strike 
entry  ana  exit  points  and  plctorlally 
represented  the  normalized  1 MHz  component 
of  the  direct  strike  akin  current  density 
distribution.  Vectors  are  drawn  at  the 
centroid  of  each  visible  triangular  face  to 

indicate  the  direction  of  Re{jc  e^Ut} 

— O 

for  a particular  cot,  where  the  skin  current 

phasor,  Jg,  is  found  at  each  centroid  from 

the  coefficients  of  triangular  patch 
expension  at  adjacent  edges.  For  the 
pictures  shown  here,  cot  was  taken  to  be 
zero.  Current  magnitude  information  may  be 
obtained  from  Fig.  5 and  7 where  numbers 
from  0 to  9 are  printed  on  each  triangular 
face  to  indicate  the  magnitude  of  the 
current  at  that  face  scaled  relative  to  the 
maximum  current. 

At  level  B of  the  interaction  table  the 
skin  currents  Jg  and  charges  Pg,  which 

are  obtained  from  Jg  by  the  equation  of 

continuity,  drive  a model  for  the  particular 
POE  under  consideration.  This  may  be  a 
composite  avionics  bay  door,  or  a joint  or 


seam  in  the  aircraft  skin.  For  each  case, 
we  relate  internal  sources  J*nt,  M*nt 
to  Jg  by: 

® Joints  - joint  admittance  Y^ 


..int 

pi 


- Y^J. 

j -» 


(3) 


° Composite  Panels  - transfer  impedance  Zt 

Mlnt  “ Zt  is  (*> 

TpQE(f)  thus  defined  by  or  Y^1 

depending  on  the  POE  under  consideration. 
The  effects  of  more  than  one  POE  are 
linearly  superposed. 

The  interior  sources  radiate  in  the 
presence  of  an  interior  compartment  to 

produce  the  interior  fields  J5'*'nt  and  H*nt. 

V I 

The  transfer  functions  T_  and  T„ 

account  for  the  effect  of  the 
compartment  detail.  An  example  of  the 
effects  of  the  compartment  transfer  func- 
tions is  shown  in  Fig.  8 where  a two- 
dimensional  field  mapping  of  the  F-14  cock- 
pit is  illustrated  for  a nose-tall  strike. 

At  successive  levels  of  interaction, 
the  complexity  cf  the  problem  car. 
exponentially  Increase  if  one  tries  to 
accurately  model  the  excitation  ot 
convoluted  masses  of  cable  bundles.  Here 
we  assume  that  a bundle  of  wires  routed 
from  one  box  to  another  within  the 
penetrated  compartment  of  interest  may  be 
modeled  by  a transmission  line  where  a 
single  TEM  mode  is  assumed.  The  interior 
fields  are  used  to  determine  distributed 
voltage  and  current  excitations  along  the 
line  so  that  we  may  write 


G^f.l)  l4<f,t)  + cJV.t) 
t)  I <f,t>  ♦ cjv(f,t) 


(5) 


The  functions  G are  obtained  from  a 

B 

straightforward  solution  of  the  trans- 
mission line  equations  using  impulsive 
voltage  and  current  sources. 

In  some  cases,  the  bulk  voltages  and 
currents  may  be  used  to  define  an  avionics 
box  protection  function  Tp  giver,  by 


(6) 
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where  I_  and  V_  are  the  threshold  box 
B B 

current  and  voltage  damage  levels, 
respectively.  Tp  gives  the  additional  box 

protection  required  against  the  computed 
bulk  voltage  V_  and  current  I_. 

B B 

CONCLUSION 

The  utility  of  the  above  procedure  is 
made  apparent  when  put  in  the  form  of  a 
user-oriented  design  program.  Data  bases 
which  describe  the  various  frequency  domain 
transfer  functions  are  accessed  as  needed 
and  can  be  updated  or  changed  depending 
upon  design  constraints.  One  e the  exact 

location  of  a cable  bundle  is  described 
inside  the  aircraft,  as  well  as  all  of  its 
associated  POEs,  induced  current  and 
voltage  levels  may  be  computed  for  a given 
lightning  stroke  situation.  An  example  of 
this  type  of  computation  for  a "generic" 
cable  harness  in  the  nose  wheel  well  of  an 
F-14  aircraft  is  given  in  Fig.  9. 
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Equipotential  lines  at  cross  section  of  F-14 
fuselage  station  320  inside  aircraft  test 
fixture.  The  fuselage  carries  the  injected 
current,  which  is  returned  to  the  generator 
by  return  conductors.  This  is  used  to 
estimate  field  levels  inside  large  apertures 
and  the  total  inductance  of  test  fixtures. 


x'8  indicate  aircraft 
test  fixture  return 
conductors 

m 83.39V  (potential 
of  fuselage) 

$r  “ “16.61V  (potential 
of  return  conductor 


A 5 -10V 


S = 80V 


5V  Intervals 


| 


i 


I 


l 


I 

i 


k 


I 

i 


I 

t 

l 

I 


I 


1 


4 

“ I 


i 


Fig.  8 - Field  Lines  Inside  F-14  Cockpit  for  Nose-Tail  Strike 
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Example  of  transmission  line  model  applied 
to  F-14  wire  bundle  coupling  in  large 
aperture  region.  Distributed  POE  fields  are 
computed  by  two-dimensional  field 
penetration  results. 
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Open-circuit  voltage  at  terminal  1 of  Bundle 
2 when  Terminal  2 is  terminated  in  50ft. 
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Fig.  9 - Generic  Cable  Bundle  Results  for  F-14 
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